Critical behavior of the 3D-Ising model on a poissonian random lattice 



in 
o 
o 

(N 
<D 
(N 



S3 

S3 
i 

C/3 



I 

S3 

o 
o 



> 
m 
o 
m 

OS 

o 
in 
o 

03 



i 

S3 
O 

o 



X 



F. W. S. Lima 1 , U. M. S. Costa 2 , and R. N. Costa Filho 2 
1 Departamento de Fisica, Universidade Federal do Piaui, 57072-970 Teresina, Piaui, Brazil 
2 Departarnento de Fisica, Universidade Federal do Ceard, 60451-970 Fortaleza, Ceard, Brazil. 

The single-cluster Monte Carlo algorithm and the reweighting technique are used to simulate the 
3D-ferromagnetic Ising model on three dimensional Voronoi-Delaunay lattices. It is assumed that 
the coupling factor J varies with the distance r between the first neighbors as J(r) oc e~ ar , with 
a > 0. The critical exponents j/v, [3/v , and v are calculated, and according to the present estimates 
for the critical exponents, we argue that this random system belongs to the same universality class 
of the pure three-dimensional ferromegnetic Ising model. 
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I. INTRODUCTION 

The Lenz-Ising model has a vast number of applica- 
tions ranging from solid-state physics to biology. It is 
the oldest and most simple model for cooperative behav- 
ior which shows spontaneous symmetry breaking^, 0- 
The symmetry breaking is very important feature of tran- 
sitions that can take place between phases with differ- 
ent symmetry like solid- liquid and ferromagnetic systems. 
The ferromagnetic transition is an example of continuous 
phase transition where the spin orientation symmetry is 
broken in the ferromagnetic phase due to the formation of 
magnetic domains. Inside each domain, there is a mag- 
netic field pointing in a fixed direction chosen sponta- 
neously during the phase transition. 

One of the main characteristics of a continuous transi- 
tion is the existence of critical exponents in the behavior 
of physical properties near the critical temperature. For 
example, the heat capacity obeys a power law 



C~ \T C -T\ 



(1) 



where a is the critical exponent associated with the heat 
capacity. When < a < 1, the heat capacity diverges 
at the transition temperature. This is the behavior in 
the 3-dimensional Ising model for ferromagnetic phase 
transition. A natural question that follows is whether 
the critical exponent is modified in the presence of im- 
purities or a random media. In such systems the influ- 
ence of quenched disorder can be classified by the Harris 
criterionj^, |j, |fj.|G( which predicts that a perturbation is 
relevant if a > 0, where a is the critical exponent of the 
specific heat of the ideal system. For a < quenched 
disorder is not relevant, while no prediction can be made 
in the marginal case a = 0. 

Several studies using Monte Carlo simulations have 
shown that for site-dilution and bond-dilution disor- 
der is relevant for the 3D ferromagnetic Ising model. 
However, when considering connectivity disorder using 
a Voronoi-Delauney lattice in the Ising 3D model, Jankc 
et.al. have shown that this kind of disorder does 
not modify the critical exponents. This result contradicts 
the Harris criteria. The explanation is that the disorder 
in the lattice was weak and more time running simula- 



tions was necessary to obtain some meaningful shift from 
the critical exponents on a regular lattice. 

In this piece of work we investigate a 3D ferromagnetic 
Ising model on a Voronoi lattice using the same algorithm 
as in Ref. 7. But here, besides the fact that the coordi- 
nation number is random, we add an additional factor to 
the randomness of the Voronoi lattice. The coupling be- 
tween the neighbors spins now depends on their relative 
distance according to the expression 



Ji 



(2) 



where J is a constant and e a > enters the model as a 
parameter When a — we have the same results as 
in Ref.0, i.e., the randomness in the number of coordi- 
nation does not change the critical exponents. We show 
here that the addition of the new variable J(r) does not 
change the universality class of the Ising model. 



II. MODEL AND SIMULATIONS 

We consider the ferromagnetic Ising model in a three 
dimensional Delauney lattice where the coordination 
number varies locally and the coupling factor depends 
on the relative distance between the first neighbors ac- 
cording to Eq[3 In order to build the lattice, we first 
distribute N sites x% randomly in a cube with edge equal 
to unit. To each site x\ a polyhedric Voronoi cell Cj is 
associated. This cell contains all the points x ^ Xi whose 
Euclidian distance to d(x, Xi) is inferior to the distance 
d(x, Xi+j). A dual lattice is then constructed joining the 
Xi of each cell that has a common edge. Using such a 
construction, the lattice connectivity is random and has 
an average value equal to 15. 

The Hamiltonian to the ferromagnetic Ising model can 
be described as 



-KH = J2 J > 



ijSiSj 



(3) 



<ij> 



where K = 1/KbT, T is the temperature, Kb is the 
Boltzmann constant, and the coupling factor is a func- 
tion defined by Eq|21 The spin Si can have two values 




In L 

FIG. 1: Log-Log plot of the average magnetization < \m\ > 

versus the system size L for three different values of the cou- FIG. 2: The logarithm of the maximum susceptibility against 
pling parameter a: a = 0.0 (circles), a = 0.5 (squares), and the system size L for a = 0.0, a = 0.5, and a = 1.0. 
a = 1.0 (diamonds). 



Si = 1 or Si = —1, and the sum is evaluated in the 
first neighbors only. For simplicity, the system size used 
is defined in terms of a regular cubic lattice with side 
L = N 1 ' 3 . In our simulations we used several system 
sizes; N = 15 3 , 20 3 , 25 3 , 30 3 , 40 3 , and 50 3 . Depending on 
the system size, we have generated a different number 
of system copies: for L = 40 and 50 we have done 25 
copies; 30 copies for L = 30, and 50 copies for L = 15, 
20, and 25. To perform the simulations we used the Wolf 
algorithm^]] close to the critical temperature K c . 

We start the simulations assigning to each site of the 
Voronoi lattice a random spin value, i.e. the initial state 
is totally random, and letting the system relax (reach 
equilibrium) for 120000 to 240000 timesteps or clusters. 
After that we generate from 2.4xl0 6 to 4.8xl0 6 clusters 
(depending on the system size) to evaluate the averages. 
For each 12 clusters the energy per spin e = j| and 
the magnetization per spin m = ^ were measured 
and registered on temporary files. In the end, we ob- 
tained from 2.0xl0 5 to 4.0xl0 5 values of the energy and 
magnetization which are is analyzed using the histogram 
techniaue|l2j|. In order to calculate the thermodynamic 
properties like specific heat C^(K) — K 2 N(< e 2 > 
— < e > 2 ), for each copy classified by the index i we used 
the reweighting technique to calculate the copies average 
indicated by the brackets in the expression below, 

C(K) = [C®(K)] = ± i 1 £C®(K) (4) 

i 

The variable R means the number of copies in the simu- 
lations. The copies average is calculated over the the C (i) 
and not over the energy moments because the quenched 
averages must be calculated at free energy level instead 
of the partition function level. Finally, we determine the 



maximum value C max = C(K max ) for each lattice size 
and analyze the scale behavior of C ma x e Kc max with the 
system size. Other thermodynamic quantities like sus- 
ceptibility x®(K) = KN(< m 2 > - < \m\ > 2 ), mag- 
netic cumulants U^l (K) = 1 - (< m 2p > /3 < \m\ p > 2 
) (p = 1,2), and several derivatives related to K, such 
as d w < \m\ > /dK, and d^ln < \m\ p > /dK, can be 
calculated in the same way. Once the dependence on the 
point Kq is well known for each run through the reweight- 
ing, we can easily calculate the quenched averages of the 
quantities. 



III. RESULTS 

To estimate the critical temperature we calculated the 
magnetic cumulants Uip using the histogram technique 
in combination with the finite size scaling analysis for 
three different values of the parameter a that controls 
the strength of the coupling between the spins. For lat- 
tices sizes ol N > 8.000 we get K c » 0.0724128 and 
C/ 2 * w 0.58567, K c » 0.0724215 and C/ 4 * w 0.4645536 
for a = 0.0, K c » 0.139265 and C/ 2 * « 0.58964, K c « 
0.139254 and [7 4 * w 0.46976 for a = 0.5, K c w 0.254941 
and C/ 2 * w 0.587601, K c w 0.254953 and U\ « 0.46765 
for a = 1.0. These results show that the universal Binder 
cumulant t/^ is different from t/j* for each value of the pa- 
rameter a but they do not change for different values of a, 
and the critical temperature is the same for both values 
of the magnetic cumulant for each a. Also, the criti- 
cal temperature decreases as the parameter a increases. 
Such behavior is very similar to the one in dilution of 
pure systems. 

To estimate the ratios of the exponents f3/v e 7/f, we 
use the fact that the magnetization on the inflexion point, 
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FIG. 3: Log-log plot of the maximum magnetization deriva- 
tive against the system size L for three different values of a: 
a — O.O(circles), a = 0.5 (squares), and a = 1.0 (diamonds). 
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FIG. 4: Log-log plot of the maximum average squared magne- 
tization derivative (EqEJ against the system size L for three 
different values of a: a = O.O(circles), a — 0.5 (squares), and 
a = 1.0 (diamonds). 



and the susceptibility on K*- nax scales with 

| < \m\ > \ lnf = L-P/vfitL 1 '") oc L-M" (5) 



XmaAL) = x(K^ ax (L),L) - ALT'/", 



(6) 



for large enough system sizes. Moreover, since the deriva- 
tives of the magnetization with respect to K depends on 
the factor IP' V which is an argument of the scale function 

/El, 



d<\m\> 
dK 1 



= L-^+V'f {tL l ' v ) oc L^ v , (7) 



the maxima of the logarithmic derivative scales like 
! din < \m\ > ! ! d < \m\ > /dK : l , v 

m.n.x OC Ij , 



dK 



and 



din < to > 



dK 



< \m\ > 
d < m 2 > /dK 



(8) 



< TO Z > 



otL 1,u . (9) 



The above expression determines the value of the expo- 
nent v. Once finding that exponent, it is possible to 
calculate the values a/v of the specific heat. Through 
the finite size scaling analysis, it is possible to argue that 
the specific heat behaves like 

C max {L) = C[Kg ax {l),l] = B + B x lnL (10) 

in the second order phase transitions. To estimate the 
exponents ratio a/v we use the hiperscale relation a/v — 
2/v-D. 

In Figure (1) we show in a logarithmic scale a graph of 
the magnetization against the lattice size L for three dif- 
ferent values of the coupling parameter a = 0.0, a = 0.5, 



and a = 1.0. The linear adjustment of the square minima 
of the simulation data gives us (3/v = 0.48920 ± 0.00823, 
13/ v = 0.51628 ± 0.00895, and f3/v = 0.51455 ± 0.00683, 
respectively. In Fig. (2) we plot the logarithm of the 
maximum susceptibility against the system size for a = 
O.O(circles), a — 0.5 (squares), and a = 1.0 (dia- 
monds). This curve inclination defines the exponents ra- 
tio j/v = 1.9223 ± 0.0295, 7/f = 1.9537 ± 0.0043 and 
7/V = 1.9758 ± 0.0069, respectively. Figs.(3) and (4) 
show a log-log plot of the maximum derivatives given by 
Eqs. IBland El respectively, against the lattice size L for 
a = 0.0, a = 0.5 c a = 1.0. The linear fitting inclination 
gives l/v 1 = 1.5184 ± 0.0302, 1/z/j = 1.5581 ± 0.0265 
and = 1.5364 ±0.0095 and l/v 2 = 1.5207 ± 0.0297, 
l/v 2 = 1.5649 ± 0.0332 and l/v 2 = 1.5326 ± 0.0078, for 
the same values of a as in the previous cases. In Fig. 5 
we show the maximum value of the specific heat against 
the lattice size L and through the hyperscale relation we 
get a/v « 0.0391, 0.123, and 0.069, for a = 0.0, 0.5, and 
1.0 respectively. 



IV. CONCLUSION 

We have calculated here the critical exponents 7/^, 
(3/v e v to the ferromagnetic Ising model in a three di- 
mensional random Delauney lattice with a coupling fac- 
tor depending exponentially on the distances between 
first-neighbors sites. 

Dilution in a regular lattice was responsible for chang- 
ing the universality class of this model, while random 
connectivity did not affect the exponents at all. The 
introduction of some kind of dilution through a distance 
dependent coupling factor in addition to the random con- 
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FIG. 5: The maximum of the specific heat against the lattice 
size for different values of the parameter a as in the previous 
figures. 



nectivity of the lattice could increase the strength of the 
quench randomness and change the exponents. Although 
the critical temperature for the system decreases as the 
coupling factor a increases indicating a dilution effect as 
in the case of a pure lattice, the results obtained here 
point out that these two ingredients are not enough to 
drive this system away from the universality class of the 
pure 3D ferromagnetic Ising model. That can be ex- 
plained by the fact that instead of increase the random- 
ness, the kind of dilution included here is responsible for 
the weakening of the long connections. This causes the 
lattice to become more regular and therefore the model 
is not affected by the effective randomness introduced. 
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